

AN APPROACH TO THE PARAMETRIC DESIGN OF ION THRUSTERS * 


Paul J. Wilbur** 

Colorado State University 
Fort Collins, Colorado 80523, U.S.A. 


and 

John R. Beattie and Jay Hyman, Jr.** 
Hughes Research Laboratories 
Malibu, California 90265, U.S.A. 


V J 



88-080 



£ 


ABSTRACT 

A methodology that can be used to determine 
which of several physical constraints can limit ion 
thruster power and thrust, under various design and 
operating conditions, is presented. The methodology 
is exercised to demonstrate typical limitations 
imposed by grid system span- to- gap ratio, intragrid 
electric field, discharge chamber power per unit 
beam area, screen grid lifetime and accelerator grid 
lifetime constraints. Limitations on power and 
thrust for a thruster defined by typical discharge 
chamber and grid system parameters when it is 
operated at maximum thrust- to -power are discussed. 

It is pointed out that other operational objectives 
such as optimization of payload fraction or mission 
duration can be substituted readily for the thrust- 
to-power objective and that the methodology can be 
used as a tool for mission analysis. 

INTRODUCTION 

Ion thrusters are very attractive propulsion 
devices not only because of their high efficiency 
and high specific impulse capabilities but also 
because they afford designers and potential users a 
great deal of design and operational flexibility. 

It is generally recognized, however, that 
operational flexibility (e.g. throttleability) 
carries with it a price of a greater system 
complexity that may not be justifiable for many 
missions and it shouLd be noted that design 
flexibility also has its unattractive aspects. It 
can for example conjure up impressions, particularly 
in the minds of those who do not work directly with 
ion thrusters, that they are too complex. In order 
to overcome this concern, the work described herein 
has been undertaken and a framework within which ion 
thruster performance and design information can be 
presented in a simple, easily understood way has 
been sought. It is hoped that this methodology can 
be used to facilitate communication, teaching, and 
an identification of the most productive areas for 
research and development in support of a variety of 
propulsion system objectives. 

The development of this methodology might be 
pursued In terms of empirical relationships 
following the^ techniques used by Byers and Rawlin 
and by Byers but the desire to use it for teaching 
purposes coupled with the fact that basic models 
describing ion thrusters have become available since 
their work was published, has prompted us Instead to 
base the development on basic physical models. 

Where parameters are required in this development to 
accommodate non- ideal behavior, every attempt has 
been made to use parameters that are physically 
meaningful, commonly used and theoretically based. 

The objective of this paper is to present a 
methodology that can be used to establish a link 
between thruster design, operational and mission 


parameters in order to determine which of many 
possible performance -limiting phenomena will 
dominate in a particular situation. A second 
objective is W* to establish a technique that can be 
used to present information on this subject In a 
manner that ts both easy to visualize and 
understand. 

THEORY 

The cannonlcal design variables associated with 
an ion thruster are the specific impulse at which it 
operates (I )„ its beam area (A.) and its beam 
power (P_) . ''Generally, one would want to prescribe 
the specific Impulse and beam area of a device and 
then determine the beam power at which it could be 
operated under a given set of design, operational 
and mission constraints, so in this development 
specific impulse and beam area will be treated as 
Independent variables and beam power will be treated 
as the dependent variable. 

Examples of physical constraints which have 
been found to limit the beam power at which a 
thruster can operate include: 

Grid system span- to- gap ratio. This constraint 
is determined by ones ability to hold grids 
close together In an environment of thermally 
Induced distortion and in some situations 
significant electrostatic attraction. It Is 
influenced greatly by mechanical design and 
fabrication considerations and while a value of 
600 might be considered a reasonable limit for 
conventional circular grids, the use of 
lntragrld supports or non-circular (e.g. annu- 
lar or rectangular grids) could facilitate 
substantial increases above 600. 

lntragrld electric field. Excessively high 
electric fields between the screen and 
accelerator grids of a thruster result in 
electrical breakdown and an inability to 
extract an ion beam. This limit is influenced 
by such factors as the surface finish and uni- 
formity of the lntragrld spacing and typical 
values in operating thrusters have generally 
been about 2 kV/mm ^although much higher values 
have been reported. 

Discharge power per unit beam area. This is 
actually a heat transfer limit that is Imposed 
because components such as magnets, anodes or 
grids can overheat if the heat removal rate is 
Inadequate. This constraint could be formu- 
lated in terms of specific heat transfer 
limitations for particular components but for 
the illustrative purposes of this paper, it 
will be assumed that the allowable discharge 
power scales directly with beam area and that a 
limit in the range 15 to 30 kW/m 2 is 
reasonable . 
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Screen grid lifetime . The lifetime of the 
screen grid as well as other components exposed 
to the discharge plasma is limited by the 
process of ion- induced sputter erosion. 

Typically grids and other components have been 
designed to have lifetimes in the range of 10 4 
to 2 x 10 4 hours. These lifetimes are 
influenced by the materials and ions involved 
and by the discharge voltage. 

Accelerator grid lifetime. The accelerator 
grid is exposed to small currents of high 
energy, charge -exchange ions that limit its 
lifetime through sputter erosion. This grid is 
typically designed to have a lifetime in the 
same 10 4 to 2 x 10 4 hr range that the screen 
grid has. 

The above list of constraints on ion thruster design 
is not exhaustive, but they represent constraints 
that have been encountered and they can be used for 
purposes of illustration. 

A trade-off exists between the propellant and 
power requirements for an ion thruster and for this 
reason a preferred operating point or operational 
objective exists. The one selected influences the 
severity of each of the design constraints identi- 
fied above. Once an objective that requires 
operation, for example, at the point where the 
thruster will produce maximum thrust per unit input 
power or maximum payload fraction on a prescribed 
mission has been identified, one can define the 
extent of the power limits cited in the preceding 
paragraphs. For the purposes of this study it will 
be assumed that the thruster is to be operated at 
the point of maximum thrust- to -power , but it should 
be recognized that this additional input, which 
could be much more complex, might come for example 
from an algorithm that optimizes a mission 
objective . 

Mathematical Development 

The beam power (P_) produced by an ion thruster 
is expressed most simply as the product of beam 
current (J n ) and beam (or net accelerating) voltage 

tV- 

P B “ J B V B (1) 

The beam current in this equation is related to the 
peak current density being extracted through the 
grids (j ) and the beam area (Ag) through the beam 
£latness m parameter (IF) which is defined as the 
ratio of average- to -peak ion beam current density. 

s V (2) 

It Is noted that the flatness parameter can be 
calculated for a given discharge chamber using^the 
finite element technique developed by Arakawa . 


the space-charge limitations that develop there. 

The final factor in Eq. 3, the effective ion 
acceleration length (i ) , is the one that is limited 
by span- to -gap and electric field considerations and 
the imposition of these constraints will be consid- 
ered in detail next. It will be assumed that the 
value of I is the same for each aperture set over 
the entire e grid surface for this development. 

The total accelerating voltage that appears in 
Eq. 3 can eliminated in favor of the beam voltage by 
introducing the net-to- total accelerating voltage 
ratio (R) 


R - V B /V T . (4) 

The beam voltage can be related to the specific 
impulse (I ) by recognizing that the specific 
impulse is Sp defined as the thrust (F) per unit 
weight flowrate of propellant (m g 0<> ) 

- ** «,o ' (5) 

that the thrust is given by the momentum equation 

F - [■ <|J [u F t aj , ( 6 ) 

and that the ion exhaust velocity (U) is related to 
the beam voltage through the conservation of energy 
expression 

V B e - U 2 /2 . O) 

In Eq. 6 the product of propellant mass 
flowrate (m) and propellant utilization efficiency 
(^ ), which Is the first bracketed term, represents 
the flowrate of thrust producing (high velocity) 
propellant and the second bracketed terra represents 
the effective jet velocity of the beam ions along 
the thruster axis. This second term includes two 
thrust correction factors, one (F c ) that reflects 
the fact that many Ions will emegge from the grid 
system on divergent trajectories and a second (at) 
that accounts for the fact that Eq. 7 describes the 
velocity of singly charged ions only when in fact 
some multiply charged gions will generally be 
produced and extracted. 


Combining Eqs. 5, 6 and 7 the desired 
expression for beam voltage in terms of specific 
impulse is obtained 


V 


B 



2e 


la 


( 8 ) 


and this may be combined with Eqs. 1 through 4 to 
obtain 


*o r ts fell 1 8 98 

9 i 2 R 3 ' 2 l "u F c “ 

e v 


(9) 


The maximum current density capability of a 
grid set is determined by space -charge limitations 
which may be described approximately us£ng the one- 
dimensionally based Child- Langmuir law. 



The permittivity of free space (e^), the electron 
charge (a) and the ion mass (m^) appearing in this 
equation are known constants and the total 
accelerating voltage (V_) and the screen grid 
transparency to ions (^ ) are at the control of the/ 
designer and/or operator. The ion current density 
enhancement factor (^ ) is a factor that will be set 
equal to unity for the purposes of this paper, but 
it could take on values greater than unity to 
account for the increased current density capabili- 
ties of grids when ions approach the^screen grid 
plasma sheath at non-zero velocities <^r high energy 
electrons are injected into the sheath to mitigate 


Equation 9 defines the maximum power constraint 
associated with the ion extraction process for an 
ion thruster as a function of its beam area and the 
specific impulse at which it is to operate. Two 
physical constraints, one associated with the 
allowable span- to- gap ratio and the other with the 
maximum allowable electric field between the grids, 
actually evolve from this equation through the 
effective ion acceleration length (i ). The fact 
that propellant utilization appears In the equation 
also serves as a reminder that an operational 
objective must be defined before unique limiting 
values of beam power can be computed as a function 
of specific impulse and beam area for these two 
physical constraints. 

Span- to -Gap Ratio Constraint 

The allowable span-to-gap ratio (N) associated 
with traditional grid sets is the grid diameter- to - 
spacing ratio. In order to accommodate non- circular 
beam cross sections, however, this ratio will be 
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Fig. 1. Aperture System Geometry 


defined here using an equivalent beam diameter and 
the grid separation will be given by the expression 


i 

g 


H v* 

N 


( 10 ) 


The preferred ion acceleration length for use in 
Eq. 9 (whlcji is based on one -dimensional theory) 
would be the one that would yield the beam current 
density actually extracted from screen/accel grid 
aperture pairs in what is really a two-dimensional 
process. One value of 2 that should come close to 
doing this is illustrated e in Fig. 1 and examination 
of the geometry of this figure suggests 2 is 
related to the grid separation distance given by 
Eq. 10 through the expression 


Discharge Chamber Thermal Constraint 


If the power required to operate a discharge 
chamber becomes too great, such failures as those 
associated with magnet, screen grid or anode 
overheating could occur. Of course the allowable 
discharge power would scale with thruster size and a 
review of limiting discharge powers on thrusters 
having various diameters has suggested that it is 
probably the discharge power per unit beam area 
(Pd/A®) that defines this constraint. 

The discharge power can be expressed in terms 
of the energy cost of a beam ion («_) through the 
equation B 


Combining Eqs. 1, 



(U) 


(14) 


The beam ton energy cost appearing in. this equation 
can be computed using Brophy's model 1 in the form 


1 -exp 


K (i -v) 


f V, 


(15) 


where the baseline plasma ion energy cost (ep) can 
be computed when the propellant and discharge 
voltage (V^) are prescribed and the extracted ion 
fraction (fl) and the fraction of the ions produced 
that go to cathode potential surfaces (f ) can be 
determined when the magnetic field and electron 
source location are specified. The parameter B is 
given by ° 


B - m C 
o o 


4 o' l 
2 



(16) 


~ / (c s + V 2 + d s /4 • < ll > 

It should be noted here that Eq. 11 differs from the 
traditions^ equation for 2 used to compute 
perveances in that it accounts for Che screen grid 
thickness (t ). If the screen grid thickness is 
left out of this expression, the implication is that 
the screen hole sheath positions itself close to the 
downstream edge of the screen grid. Aston 1 has 
shown, however, that this sheath tends to position 
Itself near the upstream edge of the hole under 
normal operating conditions and Eq. 11 is therefore 
considered correct. 

Equations 9, 10 and 11 can be combined to 
define the maximum beam power at which a thruster 
with a beam area A® and a span- to -gap ratio N can be 
operated at prescribed specific impulses provided an 
operational objective, a screen grid thickness (t ) 
and a screen hole diameter (d ) are specified. For 
this illustrative study these latter two parameters 
have been defined by specifying a grid separation- 
to-screen hole diameter ratio (2 /d ) and screen 
grid thickness of unity and 0.0005 m respectively. 

Electric Field Constraint 


where the primary electron utilization factor (C ), 
the primary electron/propellant atom total inelastic 
collision cross section the primary electron 

containment length (2 ), the propellant atom 
thermal velocity (v ) and the transparency of the 
grids to neutral atoms ) may be assumed constant 
for a particular propellant and discharge voltage. 

An expression for the propellant flowrate per unit 
beam area needed in Eq. 16 can be obtained by 
combining the definition of propellant utilization 
efficiency 


.daJt 


( 17 ) 


with Eqs. 2 through 7 to obtain 
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Equations 14, 15, 16 and 18 when combined with an 
operational objective yield the beam power/beam 
area/specific impulse surface that defines the 
discharge power per unit beam area constraint. 


In order to prevent electrical breakdown 
between grids it is presumed that a limiting 
electric field E cannot be exceeded. Assuming a 
uniform grid spacing, this may be expressed 
mathematically using the expression 


2 


g 


Ii 

E 


R E 


( 12 ) 


Combining this with Eqs. 9 and 11, an operational 
objective, a screen grid thickness and a screen hole 
diameter, the beam power/beam area/specific impulse 
surface limited by electrical breakdown 
considerations is defined. 


Screen Grid Lifetime Constraint 

Although the screen grid is assumed to be the 
life limiting component subjected to sputter erosion 
in this paper, additional constraints pertaining to 
other components that are subjected to ion 
bombardment could impose similar constraints. They 
would be Incorporated into the analysis In the same 
way as the one developed here for the screen grid. 

If the screen grid lifetime is considered to have 
expired when a fraction 7 of its initial thickness 
t g has beengsputtered away then this Lifetime will 
be given by 


469 



Accel Grid Lifetime Constraint 


r 

s 


8 ' N * h 


K S X 


,) + °.s j 4 jC (2 V * 


(19) 


where p and M are the density and molecular weight 
of the grid material, N is Avagadro's number, j + 
and J are the singly and doubly charged ion 
current densities striking the grid and S^(V_) 

S (2V n ) are the sputtering yields of the grid 
material for singly and doubly charged ions 
evaluated respectively for ions with energies equal 
to the discharge voltage and twice the discharge 
voltage. Recognizing that the grid will sputter 
most rapidly at the point of maximum current density 
m - \ ) Eqs. 1, 2 and 8 can be combined with 

Eqt 19 cS* obtain 


P * 7, 


P B- 
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( 20 ) 


The doubly- to -singly charged ion current density 
ratio (J /j ) appearing in this equation is a 
strong functlon^f discharge operating conditions 
and is given by 
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( 21 ) 


Some simplification of this equation is generally 
possible because the first term is negligible at 
typical discharge chamber electron temperat^s* and 
energies. In Eq. 21 v and v are the Bohm and 
neutral atom thermal velocities, and are the 

transparencies of the grids to ions and neutral 
atoms and Q , P and etc. are the rate factors 
associated with^production of the various ionic 
species. The primary- to- Kaxwe Ilian electron 
density ratio (n /n^) appearing in Eq. 21 is given 
by 11 P 


n 


_a 
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( 22 ) 


In this equation, the primary electron velocity (v ) 
is determined by its energy which is assumed in turn 
to be equal to the discharge voltage. The discharge 
chamber volume- to -beam area ratio (V/A_) I s simply 
the discharge chamber Length which may oe assumed to 
be relatively insensitive to beam area (e.g. V/A^ - 
0.1 + 0.1 [ 1-exp ( -lOAg) ] has been used for the 
examples presented here) . The final expression 
needed to define the screen grid lifetime constraint 
is one for the propellant mass flowrate per unit 
beam area. It is obtained by combining Eqs. 1, 8 
and 17 . 




2P n 


I g 

sp °eo 


(23) 


Equations 20 through 23 can be solved for the beam 
power/beam area/specific impulse surface that 
defines the screen grid lifetime constraint when an 
operational constraint is specified and design 
definitions associated with the chamber (e.g., pro- 
pellant, discharge voltage, etc.) have been made. 


For a three-grid ion optics system where 
sputtering on the barrel (interior surface) region 
of the accel grid holes dominates the erosion of the 
grid, the change in accel hole diameter (d^) per 
unit time is given by 


d(d a ) 

dr 


S j 


e N p w d 


7 V f M 

as a. 


(24) 


In this equation S is the sputter yield of the 
accelerator grid material at the prevailing charge 
exchange ion kinetic energy, n Q is the neutral atom 
density in the charge exchange reaction region, 
is the volume of this region, f is a factor that 
describes the extent to which these ions are focused 
or distributed along the barrel region and t^ is the 
accelerator grid thickness. For a two-grid optics 
set where sputtering on the downstream surfage of 
the accel grid dominates, a similar equation is 
used but for the example being considered here, the 
three- grid equation (Eq. 24) will be used 
exclusively. 


The neutral atom density in the charge exchange 
region is given approximately by the expression 


A B *o V o e 


(25) 


If the allowable change in diameter is Ad and is 
the desired lifetime then the integrated form of * 
Eq. 24 may be combined with Eqs. 1, 2, 8 and 25 to 
obtain 
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(26) 


In order to compute beam power as a function of beam 
area and specific Impulse using this equation the 
volume of the charge exchange region has been 
assumed to be equal to that of the cylinder shown in 
Fig. 1 with a diameter equal to that of the accel 
grid aperture (d ) and a length twice the grid 
separation distance <i ). In addition the ratio of 
accel grid thickness to®grid separation distance 
(t /i ) has been assumed Co be constant. 

Operational Objective 


It is necessary to define an operational 
objective for a thruster in addition to the various 
physical constraints imposed on it in order to 
define the thruster discharge chamber operating 
point, i.e., the propellant utilization efficiency 
and discharge power at which the thruster should be 
operated. This can be done for example by defining 
a mission of interest, which might be characterized 
by a mission time and characteristic velocity or a 
complex mission algorithm; computing the masses of 
the system elements (power plant, propellant, 
payload, etc.) and determining the propellant 
utilization efficiency operating point at which the 
associated payload fraction is maximized. The 
optimum propellant utilization operating point has 
been selected in this study so the thruster 
operates at the maximum thrust- to -power condition. 
This condition is defined by combining Eqs. 1, 5, 6 
7, and 17 with the expression for thruster 
electrical efficiency which Is given by 







Z? 


(27) 


In this equation P_ is the totaL thruster power and 
EP is the sum of tne powers needed to generate ions 
and to sustain chruster temperatures, propellant 
flow- races and neutralizer operation (i.e., the 
power loss cerm) . Frequently the dominant power 
loss is the discharge power required to produce Che 
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Iona and it can be expressed in terms of the energy 
cost of a beam Ion (P_ - «_). Assuming the 

discharge power does u dominate the losses, Eq. 27 
can be rewritten as 


and this equation can be combined with those 
identified above to obtain the expression for 
thrust -to- total power. 


Table X (continued) . 


L i 1 


I SB S SQ 

h» F t 

2 + :*i 
■i j 


£_ _ Sp g ?9 


By seeking the propellant utilization that maximizes 
this equation at each specific impulse, operation at 
maximum thrust- to-power is realized for each 
thruster beam area and physical constraint 
condition. 


The analysis technique outlined in the 
proceeding section can be used to Investigate the 
effects of a wide variety of design and operational 
parameters on the power and thrust capabilities of 
ion thrusters. Because the purpose of this paper is 
to demonstrate the capability of the methodology 
involved rather than to draw conclusions based on an 
exhaustive study conducted using it, one set of 
typical values of thruster parameters has been 
selected for use in the analysis. The values used 
are listed in Table I and while they are considered 
to be typical of ion thrusters in general, they do 
not represent any particular thruster. The rate 
factors used in the analysis are based on an assumed 
Maxwellian electron temperature of 5 eV and a 
primary electron energy, which is consistent with 
the discharge voltage (30 eV) . 

Table I. Thruster Parameters Used in Example Study 
in the Order of Their Appearance 


Parameter/Property 

Symbol 

Value Used 

Discharge Voltage 

V D 

30 V 

Total Electron- Atom Inelastic 
Collision Cross Section 16 16 

a * 
o 

, , in -20 2 
7.4x10 m 

Primary Electron Containment 
Length 

I* 

3.5 m 

Neutral Atom Thermal Velocity 

V 

0 

290 m/sec 

Allowable Screen Grid Erosion 
Fraction 

7 s 

0.5 

Ionization Rate Factors for 
5 eV Temperature Maxwellian 
Electrons 13 

< 

c 

7. lxlO'^mVsec 
5.1x10 ^m*Vsec 

Ionization Rate Factors for 
30 eV Energy Primary 
Electrons 13 

p + 

0 

c 

1.3x10 ^m 3 /sec 
2.5x10 ^nV s ec 

Sputter Yields for Molybdenum 
Screen Grid for Singly and 
Doubly Charged Ions 17 

s + 

s 

S++ 

s 

2x10* 6 
1.5xl0* 3 

Sputter Yield for Accel Grid 17 

S 

a 

1.0 

Charge Exchange Cross 
Section 13 

a 

ce 

3x10 * 19 m 2 

Accel Grid Thickness-to-Grid 
Separation Ratio 

Allowable Accel Crid Erosion 

t /t 

a g 

Ad a /d 

0.3 

0.5 

a 

Charge Exchange Ion Focusing 
Factor 

f 

a 

0.25 


Parameter/Property 


Ion Beam Flatness 

Ion Current Density 
Enhancement Factor 

Screen Crid Transparency 
to Ions 

Crid System Transparency 
to Neutral Atoms 

Ion Mass (Xenon) 

Net-to-Total Accelerating 
Voltage Ratio 

Beamlet Divergence Thrust 
Factor 

Multiply Charged Ion Thrust 
Factor 

Screen Grid Thickness 

Grid Separation to Screen 
Hole Diameter Ratio 

Baseline Plasma Ion 
Energy Cost 

Extracted Ion Fraction 

Fraction of Plasma Ions Going 
to Cathode Potential Surfaces 


Value Used 


2.2xl0* 25 kg 


5x10 *^m 


50 eV/ion 
0.5 


If one prescribes a thruster with a beam area 
of 0.2 m 2 and calculates the beam power limits for 
each of the physical constraints described in the 
proceeding section with the objective of operation 
at maximum thrust- to-power , the curves of Fig. 2 are 
obtained for the values of the constraints cited in 
the figure. For each constraint, operation below 
and to the right of the associated curve assures 
operation that does not violate the constraint. 

Hence in the case of Fig. 2 the 600 span- to -gap 
constraint limits the power that can be extracted up 
to a specific impulse of -3000 sec and then the 


0.2 rn BEAU AREA 
WAX. THRUST -TO— POWER 
30 V DISCHARGE 
XENON PROPELLANT 


10 hr ACCEL — . 
LIFE CONSTRAINT 


10 hr SCREEN ■*- 
UFE CONSTRAINT 


15 kW/m* DISCHARGE^ - - 
POWER CONSTRAINT 


V 


2 WV/mm ELECTRIC 
' FIELD CONSTRAINT 


-600 SPAN -TO -GAP CONSTRAINT 


0 1 2 3 4 5 6 

SPECIFIC IMPULSE (ksec) 

Fig. 2. Typical Power Constraint Curves - 0.2 m 2 
Beam Area (Span- to -Gap and Intragrid 
Electric Field Limiting) 
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2 kV/nm electric field constraint becomes limiting. 
For the particular parameters associated with the 
other constraints, Fig- 2 indicates none of the 
other constraints (discharge power, screen grid 
life, or accel grid life) become limiting at any 
specific impulse up to 6000 sec. The fact that 
span-to-gap is limiting at low specific impulses 
while electric field becomes limiting at higher 
ones is in qualitative agreement with experimental 
observations . 

If the grids were designed so the electric 
field limit could be increased to 4 kV/mm and the 
other constraints were held fixed at the values of 
Fig. 2, then the data of Fig. 3 are generated. They 
suggest the 600 span- to -gap limit would prevail to 
-4000 sec. specific impulse and the discharge power 
per unit beam area limit of 15 kV/m 2 would become 
constraining beyond that point. Increasing the 
electric field limit has in this case allowed the 
beam power at 4000 sec I to increase from -15 kV 
(Fig. 2) to -25 kW (Fig? P 3) and at 6000 sec from 
-27 kW (Fig. 2) to -45 kU (Fig. 3). If, on the 
other hand, it were necessary to have screen and 
accel grid lifetimes of 2 x 10 4 hr then the 
accelerator grid lifetime would become limiting over 
the I range from about 3500 sec to 4500 sec as 
the daiS of Fig. 4 show. In this case the power at 
4000 sec I would be limited by accel grid lifetime 
considerations to -21 kW. 



Fig. 3. Typical Power Constraint Curves - 0.2 m 2 

Beam Area (Span- to -Gap and Discharge Power 
Limiting) 


If the beam power limits imposed by the 
constraints being considered as functions of both 
beam area and specific impulse are sought rather 
than holding beam area constant, then beam power/ 
beam area/specific impulse surfaces like those shown 
in Fig. 5 are generated. These surfaces define 
power limits associated with each constraint 
indicated and they are all plotted on the same 
scale, namely the one defined in Fig. 5a. Operation 
is permitted at any point beneath the surface 
associated with a particular constraint, but not 
above It. All of these surfaces behave as one would 
expect with the beam powers being lowest for low 
beam areas and low specific impulses except the one 
associated with screen grid lifetime. The high 
allowable power observed at low specific impulses 
for the screen life constraint seems unusual. It Is 
a consequence of operation at maximum thrust- to - 
power which implies a low propellant utilization at 
low specific impulses (Eq. 29). Low propellant 
utilizations in turn imply low doubly- to-singly 
charged ion current densities (Eq. 21) and hence 
high allowable beam powers (Eq. 20). 



SPECIFIC IMPULSE (ksec) 


Fig. 4. Typical Power Constraint Curves - 0.2 m 2 

Beam Area (Span- to-Gap , Discharge Power and 
Accel Grid Lifetime Limiting) 

When all constraints shown in Fig. 5 are 
applied and the overall limiting surface is identi- 
fied, the one shown in Fig, 6 is obtained. As the 
labels on this surface Indicate, the span-to-gap and 
electric field constraints represent the most 
restrictive limitations (span-to-gap at low specific 
impulses and electric field at higher ones) over the 
beam area range from zero to 1 m 2 . 

If one defines a mission along with powerplant 
and power processor specific mass parameters and 
then Imposes the power limits defined by all of the 
constraints associated with the data in Fig. 5, 
payload fractions can be computed as a function of 
beam area and specific impulse. Figure 7 shows the 
payload fraction surface computed when the most 
restrictive of these constraints are applied for a 
10 4 m/sec mission to be accomplished in 0.67 year 
when the powerplant and power conditioner are 
characterized by the specific masses and the power 
conditioner efficiency values cited on the figure. 
The payload fraction reaches a relatively flat peak 
at a value of 55 percent when the specific impulse 
Is near 4500 sec for all but the smallest beam 
areas. The payload fraction peak in Fig. 7 is seen 
to be very broad, so the thruster could be operated 
at specific impulses ranging from -3000 to 6000 sec 
and deliver about the same payload fraction. 

The thrust at which the ion thruster operates 
can also be computed as a function of beam area and 
specific impulse by using Eqs. 28 and 29 in con- 
junction with the most restrictive beam power 
constraining surface data (Fig. 6) . Figure 8 shows 
how the thrust associated with the constraints 
defined In Fig. 5 would vary as a function of beam 
area and specific Impulse. 

If the electric field constraint is increased 
to 4 kV/mm, the results of Fig. 3 showed that the 
discharge power per unit beam area became the power- 
limiting surface at high specific impulses when the 
beam area was 0.2 m 2 . The effect of Introducing 
this new electric field constraint for beam areas 
ranging from zero to 1.0 m 2 is shown in Fig. 9. 
Comparison of the data In Figs. 6 and 9 shows that 
increasing the electric field limit facilitates a 
substantial increase in beam power at high beam area 
and specific impulse values (250 kv vs. 150 lew) . 

If the grid lifetime constraints are tightened 
by requiring 20.000 hr operating times, then the 
accel grid lifetime constraint identified in Fig. 10 
limits the beam power in the moderate specific 
impulse- low beam area regime. Away from this region 
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Fig. 9. Typical Overall Power Constraining Surface 
Fig. 8. Typical Thrust Surface (Span-to-Gap and Discharge Power Limiting) 
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Fig. 10. Typical Overall Power Constraining Surface 
(Span- to -Gap , Discharge Power and Accel 
Grid Life Limiting) 


Fig. 11. Constant Power Contour Diagram (Span- to - 
Gap, Discharge Power and Accel Grid Life 
Limiting) 


the power limit is unchanged from the values 
presented in Fig. 9. 

The three-dimensional plots shown in Figs. 5 
through 10 are useful to illustrate qualitative 
behavior, but they are difficult to use as a source 
of quantitative information. Quantitative infor- 
mation can be presented better in the form of an 
equal beam power contour plot like the one in 
Fig. 11. This figure, which presents the same data 
as that in Fig. 10, clearly indicates the beam power 
limit at any beam area and specific impulse point. 

It also shows the constraint that is preventing 
operation at higher power levels at each beam area 
and specific impulse. Similar figures could, of 
course, also be generated to show the limiting 
values of thrust and payLoad fraction as a function 
of thruster beam area and specific impulse. 

FUTURE DIRECTIONS 

The intent of the preceding discussion has been 
1) to present a methodology and framework within 
which gross parameters describing ion thruster 
behavior could be used to predict ion thruster per- 


formance limits, 2) to cite references that describe 
how these gross parameters can be computed, and 
3) to suggest how results obtained from the analysis 
might be presented In a readily understood format. 
The analysis should, however, not be considered 
fully developed. The following statements describe 
changes that might be introduced to improve the 
analysis . 

1) The beamlet divergence (thrust) factor (F^), 
the grid separation- to- screen hole diameter ratio 

(i /d ) and the net- to- total acceleration voltage 
ratio S (R) have all been treated as constants in the 
cases presented. These quantities are, however, 
variable and they are related to each other. One 
could therefore use data like those in Ref. 19 to 
determine values of F t for prescribed values of 
i /d and R or one could incorporate additional 
operational objectives that would assure the 
parameters are selected to maximize a thruster 
parameter of Interest. 

2) The screen and accel grid thicknesses (t g 
and t respectively) have been treated as constants, 
and depending on the beam geometrical cross section 
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it may be that they would be better represented as 
functions of bean area. For example, a functional 
relationship between these thicknesses and the beam 
area based on mechanical deformation considerations 
might be desirable. 

3) Although some variation in discharge chamber 
plasma properties has been allowed in the develop- 
ment of the screen life constraint, more can be done 
in modeling the other constraints to reflect the 
effects of changes in Maxwellian electron tempera- 
ture, primary and Maxwellian election densities, the 
baseline plasma ion energy cost (e ) and the doubly- 
charged ion thrust factor (a) induced by changes in 
propellant utilization and discharge voltage. It is 
noted that the physically based models needed to do 
this are available. 

4) The grid spacing has been assumed constant 
across the entire intragrid region. The effects of 
variable spacing induced by grid thermal distortion 
and electric field induced deflection forces may 
cause this spacing to change as a function of 
location on the grids. This in turn will influence 
the ion extraction and electrical breakdown 
capabilities of the grids. 

Finally, it is noted that many of the 
parameters needed to model the plasma discharge and 
ion extraction phenomena involved in this analysis 
can be computed from basic principles. The note- 
worthy exception to $his is the primary electron 
containment factor 2 . It can be Inferred from 
discharge chamber tests, but a model that can be 
used to calculate it for a discharge chamber having 
prescribed dimensions and magnetic field character- 
istics is needed. Additional work may also be 
needed to describe more accurately the energy of 
primary electrons extracted from a hollow cathode as 
a function of discharge voltage. 

CONCLDS IONS 

A methodology that can be used to determine the 
maximum power level at which an ion thruster can be 
operated and a framework within which the resultant 
data can be presented has been developed. Physical 
constraints associated with the allowable grid span- 
to-gap ratio, the intragrid electric field, the 
discharge power per unit beam area and the screen 
and accel grid lifetimes have been identified as 
power (or thrust) limiting and relationships that 
quantify each of the constraints have been 
presented. When the methodology is exercised for a 
thruster operating on xenon at propellant utiliza- 
tion efficiencies that Induce max imum thrust- to- 
power with constraints of 600 span-to-gap ratio, 

2 kV/mm electric field, 15 kW/m 3 discharge power per 
unit beam area and 10 4 hr screen and accel grid 
lifetime requirements, the span-to-gap constraint is 
shown to be limiting at low specific impulses and 
the electric field is shown to be limiting at higher 
ones. If the electric field limit is increased to 
4 kV/mm, the discharge power becomes limiting at 
high specific impulses and if the grid lifetime 
requirement is increased to 2 x 10 4 hr then accel 
grid erosion can become limiting at intermediate 
specific impulses. Although the input data have not 
been selected so results predicted in the analysis 
can be compared to experimental observations, the 
general trends appear to be consistent with 
generally observed experimental trends. 
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